Observation of decays to p(770)7r and p(2150)7r 
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tp' — > p(770)7r is observed for the first time in a data sample of 14 million ■(/)' decays collected by 
the BESII detector at the BEPC. The branching fraction is measured to be B{ip' p(770)7r) = 
(5.1 ± 0.7 ± 0.8) X 10^^, where the first error is statistical and the second one is systematic. A 
high mass excited p state with mass around 2.15 GeV/c^ is also observed with B{ip' —* p(2150)7r 
TT+TT^Tr") = (19.4 ± 2.5t.2\^) ^ 10~^. The branching fraction of ip' tv'^tv~tv''' is measured with 
improved precision, B{ip' tv'^tv~tv^) = (18.1 ± 1.8 ± 1.9) x 10~^. The results may shed light on 
the understanding of the longstanding "pn puzzle" between J/ip and tp' hadronic decays. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 
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From perturbative QCD (pQCD), it is expected that 
both J/f/; and ■0' decaying into light hadrons are domi- 
nated by the annihilation of cc into three gluons or one 
virtual photon, with a width proportional to the square 
of the wave function at the origin Q. This yields the 
pQCD "12% rule", that is 
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■4i'~ 
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12%. 



A large violation of this rule was first observed in decays 
to pn and K*+K^ + cc. by Mark II 0, the so called 
pn puzzle. Since then BES has measured many two-body 
decay modes of the tp'; some decays obey the rule while 
others violate it 0]. There have been many theoretical 
efforts trying to solve the puzzle 0. However, none has 
been accepted as the solution to the problem. 

In the study of the pir puzzle, ip' — > pn is one of the 
key decay modes and is of great interest to both theo- 
rists and experimentalists. A recent calculation of the 
-0' — > pn branching fraction, done in the framework of 
SU(3) symmetry, takes into consideration interference 
between ^p' resonance decays and the continuum ampli- 
tude 1^; a branching fraction ofip' PT^ around 1 x 10"'* 
is predicted with a large error due to the limited precision 
for if)' decays into other vector pseudoscalar (VP) modes. 
The measurement of the ij}' ~^ P"" mode is a direct test 
of many models proposed to solve the pir puzzle 0, . 

The data used for this analysis are taken with the Bei- 
jing Spectrometer (BESII) detector at the Beijing Elec- 
tron Positron Collider (BEPC) storage ring operating 
at the ip' energy. The number of 0' events is 14 ± 0.6 
million determined from the number of inclusive 
hadrons, and the luminosity is (19.72 ± 0.86) pb~* as 
measured by large angle Bhabha events. 

BESII is a conventional solenoidal magnet detector 
that is described in detail in Refs. 0jl3- A 12-layer ver- 
tex chamber (VC) surrounding the beam pipe provides 
trigger information. A forty-layer main drift chamber 
(MDC), located radially outside the VC, provides tra- 
jectory and energy loss (dE/dx) information for charged 
tracks over 85% of the total solid angle. The momentum 
resolution is Cp/p — 0.017-\/l + p'^ {p in GeV/c), and the 
dE/dx resolution for hadron tracks is ^ 8%. An array of 
48 scintillation counters surrounding the MDC measures 
the time-of-flight (TOF) of charged tracks with a resolu- 
tion of ^ 200 ps for hadrons. Radially outside the TOF 
system is a 12 radiation length, lead-gas barrel shower 
counter (BSC). This measures the energies of electrons 
and photons over ~ 80% of the total solid angle with 
an energy resolution of as/E — 22%/ Ve {E in GeV). 
Outside of the solenoidal coil, which provides a 0.4 Tesla 
magnetic field over the tracking volume, is an iron flux 
return that is instrumented with three double layers of 
counters that identify muons of momentum greater than 
0.5 GeV/c. 



A phase space Monte Carlo sample of 2 million ■0' — * 
7r+7r~7r'^ events is generated for the efficiency deter- 
mination in the partial wave analysis (PWA). Monte 
Carlo samples of Bhabha, dimuon, and inclusive hadronic 
events generated with Lundcharm |^ are used for back- 
ground studies. The simulation of the detector uses a 
Geant3 |T3 based program, which simulates the detector 
response, including the interactions of secondary parti- 
cles with the detector material. Reasonable agreement 
between data and Monte Carlo simulation has been ob- 
served in various channels tested, including e+e^ — > 



(7)e^ 



(7)/^+/^ , J/ip ^ pp, and tp' 



J/V'TT+TT-, J/tP £+£- {i = e, 

The final state of interest includes two charged pions 
and one neutral pion which is reconstructed from two 
photons. The candidate events must satisfy the following 
selection criteria: 

1. A neutral cluster is considered to be a photon can- 
didate when the deposited energy in the BSC is 
greater than 80 MeV, the angle between the near- 
est charged track and the cluster is greater than 
16°, the first hit of the cluster is in the beginning 
six radiation lengths of the BSC, and the angle be- 
tween the cluster development direction in the BSC 
and the photon emission direction is less than 37°. 
The angle between two nearest photons is required 
to be larger than 7°. The number of photon candi- 
dates after selection is required to be two. 

2. There are two charged tracks in the MDC with net 
charge zero. A track should have a good helix fit 
and satisfy | cos 6*1 < 0.80, where 9 is the polar angle 
of the track in the MDC. 

3. For each charged track, the TOF and dE/dx mea- 
surements are used to calculate values and the 
corresponding confidence levels for the hypotheses 
that the particle is a pion, kaon, or proton (Pro6,r, 
Prohx, Probp). At least one charged track is re- 
quired to satisfy Prober > PtoBk and ProbT^ > 
Probp. Radiative Bhabha background is removed 
by requiring the tracks have small dE/dx or small 
energy deposited in the BSC. Dimuon background 
is removed using the hit information in the muon 
counter. 

4. A four-constraint kinematic fit is performed under 
the hypothesis ip' — > 7777+ 7r~, and the confidence 
level of the fit is required to be greater than 1%. 
A Four-constraint kinematic fit is also performed 
under the hypothesis of ip' '^jK^K" , and 
X^/^ttt: < X^-yKK is required to remove K+K^Tr° 
events. 

5. To remove background produced by '0' cascade de- 
cays to J/V' with J/ip — > fi^ fi^ , the invariant mass 
of 7r+7r~ is required to be less than 2.95 GeV/c^. 
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After applying the above selection criteria, the invari- 
ant mass distribution of the two photons is shown in Fig- 
ure QJi. A clear tt" signal can be seen. A fit to the mass 
spectrum (shown in Figure ^) using a 7r° signal shape 
determined from Monte Carlo simulation and a polyno- 
mial background yields 260 ± 19 7r°s. 
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FIG. 1: Two photon invariant mass distribution after final 
selection for (a) -i/)' data and (b) continuum data. The his- 
tograms are data, and the curves show the best fits. 

The contribution from the continuum llll | is mea- 
sured using (6.42 ± 0.24) pb^^ of data taken at \/s = 
3.65 GeV (continuum data). Figure shows the 77 in- 
variant mass distribution and the fit. The number of tt*^ 
from the fit (10.0 ± 4.2) is subtracted incoherently from 
the 1}}' data after normalizing by the two luminosities. 
This yields 229 ± 23 observed ij]' 7r+7r~7r'' events. 

Dalitz plots of the tt+tt^tt^ system for the ip' and con- 
tinuum data are shown in Figure|3 after requiring the in- 
variant mass of the two photons lies within ±30 MeV/c^ 
of the nominal tt'^ mass. (The mass resolution is around 
17.5 MeV/c^ from Monte Carlo simulation.) For the -0' 
sample, 250 events are obtained with about 13% non-7r° 
background, while for the continuum sample, 11 events 
are obtained with about 42% non-Tr*^ background. In t/j' 
decays, besides clear p bands at the edges of the Dalitz 
plot, there is a cluster of events in the center. This is 
very different than the Dalitz plot for J/tj} ^ 7r"'"7r~7r° 
decays |l2| . indicating different decay dynamics between 
J /ip and i/"' ^ tt+tt^tt''. There is no clear intermediate 
state in the continuum data. 
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FIG. 2: Dalitz plots of tt^tt tt" for (a) ip' data and (b) con- 
tinuum data after the final selection. 



The selected events are fitted in the helicity amplitude 
formalism with an unbinned maximum likelihood method 
using MINUIT For the process 

V'(l-) ^ p(l-) + 7r(0-) 

'-^ 7r(0-) + 7r(0-), 

the intensity distribution dl for the final state is written 
as 

dl= Y,{\A,\^ + \C.,\^)d{LIPS), 

i=±l 

where Ci is an incoherent background term, that is as- 
sumed to be either a constant or to have the same angu- 
lar distribution as Ai. The difference between these two 
fits is taken as the systematic error on the background 
description. LIPS denotes the Lorentz-invariant phase 
space, and the amplitude 

= ^0 (tt- , ^+ ) + (7r+ , tt" ) + A,- (vr" , ^- ) , 

where i = +1 or —1 is the helicity of the ■0', the first pion 
in each set of parentheses is the "designated" pion, and 

A^;^ — i3(m^) sin 6*^ (cos ± i cos sin (/)^)e''^*'^. 

Here c = 0,+l,or— lis the net charge of the dipion sys- 
tem, 6 and <j) are the polar and azimuthal angles of the 
p in the ip' rest frame, Ot^ and ipT, are the polar and az- 
imuthal angles of the designated pion in the p rest frame, 
and B{m?) describes the dependence of the amplitude on 
the dipion mass m: 

""^^ ^ ^ ttt;^^ ' 

where, BW{m^) is the Breit-Wigner form of the p(770) 
or its excited states. Here, the Gounaris-Sakurai param- 
eterization is used; j3j and Cj are the relative phase 
and the relative strength, respectively, between the ex- 
cited p state j and the p{770). 

Since the number of events is limited, the masses and 
the widths of all states in the fit are fixed to their PDG 
values 01 , and the number of background events is fixed 
to the number determined from the 77 invariant mass fit. 
A fit with p(770), p(1450), p(1700) and p(2150) results in 
insignificant p(1450) and p(1700) contributions. The fit 
after removing these two components yields a likelihood 
decrease of 10.7 with four less free parameters. The fit 
results are shown in Figure |31 the fit describes the data 
reasonably well. 

The fit parameters and results are given in Table 
where for results without errors, the parameter is fixed. 
The fit yields (28 ± 3)% p(770)7r in all 7r+7r-7r° events. 
By comparing the likelihood difference with and without 
the p{770)7T in the fit, the significance of yo(770)7r is deter- 
mined to be 7.8<T. The significance of p(2150)7r is larger 
than 10(7. 



4 



o 

>60 



0) ' 

> 



- p(2150) 

- P(770) 
HistogramlM 

• Data 




FIG. 3: Comparison between data (dots with error bars) 
and the final fit (solid histograms) for (a) two pion invari- 
ant mass, with a solid line for the p(770)7r, a dashed line for 
the p(2150)7r, and a hatched histogram for background; (b) 
the p polar angle in the ip' rest frame; and (c) and (d) for the 
polar and azimuthal angles of the designated tt in p helicity 
frame. 



TABLE I: t]/ Tx'^'K tx° fitting parameters and results. For 
the numbers with no errors, the values are fixed in the fit. 



Quantity 



Fit result 



Mp(77o)(GeV/c^) 
rp(770)(GeV/c2) 

/?p(2150)(°) 

Mp(2i50)(GeV/c2) 
rp{2i50)(GeV/c2) 

B(7r+7r-7r°):B(p(770)7r):B(p(2150)7r) 



0.7711 
0.1492 
-102 ± 10 
2.149 
0.363 

1:0.28 ± 0.03:1.07 ±0.09 



The fit quality is checked using Pearson's test by 
dividing the Dalitz plots into small areas with at least 
20 events and comparing the number of events between 
data and normalized Monte Carlo simulation. A x^/ndf 
= 14.6/7 = 2.1 is obtained, which corresponds to a con- 
fidence level of 4%. A fit with the p(2150) width free; or 
a fit with /9(770), /9(1450), /9(1700), and p(2150); or even 
with an extra excited p state does not improve the fit 
quality significantly. These cases show that the change 
in the number of p(770)7r events is less than 9.1%, which 
is included as part of the systematic error. The number 
of p(2150)7r events increases by 57% when other excited 
p states are added in the fit due to the interference; this 
is also included in the systematic error. 

Using the parameters of the fit in the Monte Carlo 
generator, the efficiency of if)' tt~^tt~tt° is estimated to 
be 9.02%, and the corresponding efficiencies for p(770)7r 



and p(2150)7r are 10.54% and 8.70%, respectively. 

Systematic errors in the ip' ~* 7r+7r~7r° branching frac- 
tion measurement come from the kinematic fit, the MDC 
tracking, charged particle identification, photon identifi- 
cation, background estimation, continuum subtraction, 
etc. All sources considered are listed in Tabled Most of 
the errors are measured using clean exclusive J/ip and ip' 
decay samples , while some others were described 

above. For the jo(770)7r and p(2150)7r, the uncertainties 
of fitting with different high mass p states, etc. are also 
included. The total systematic error for ijj' — ^ 7r+7r~7r'^ 
is 10.5%, and those for -0' ^ p(770)7r and /3(2150)7r are 
16.0% and respectively. 



TABLE II: Summary of systematic errors on B{tp' 
tt+tt-tt"). 



Source 


Relative error (%) 


Trigger 


0.5 


MDC tracking 


4.0 


Kinematic fit 


6.0 


Photon efficiency 


4.0 


Number of photons 


2.0 


Background estimation 


3.6 


Particle ID 


negligible 


Total number of V'' 


4.0 


Continuum subtraction 


3.0 


Total 


10.5 



Using the numbers obtained above, the branching frac- 
tions of tp' — > 7r+7r~7r'', p(770)7r and p(2150)7r are 



B(7r+7r"7r°) = (18.1 ± 1.8 ± 1.9) x 10^^ 



S(p(770)7r - 
S(p(2150)7r 



(5.1 ±0.7 ±0.8) X 10~^ 
(19.4±2.5t"i^) X 10"^ 



where the first errors are statistical and the second sys- 
tematic. 

Our B{ij' 7r+7r-7r°) agrees with the Mark II 
result within 1.8cr, and B{i/j' — > p(770)7r) is below the 
Mark II upper limit and in agreement with one model 
prediction |5j. It should be noted that the continuum 
amplitude which is considered incoherently in this anal- 
ysis could increase the p(770)7r branching fraction due 
to interference with the resonance j^. This should be 
considered in a higher statistics experiment. 

Comparing with the corresponding J/ip decay branch- 
ing fractions, it is found that both tt+tt^tt^ and p{77Q)tt 
are highly suppressed compared with the "12% rule", 
while for p(2150)7r, there is no measurement in J/ip de- 
cays. It could be enhanced in ip' decays since the phase 
space in J/ip decays is limited due to the large mass of 
the excited p state. It should be noted that using the 
J/ip and tp' — > piT branching fractions, the ip" — > pir 
branching fraction and the e"'"e~ pir cross section at 
y/s — 3.773 GeV can be predicted in the S- and D-wave 



5 



mixing model [I3 , which is proposed as a solution of the 
pvr puzzle in if)' decays. 

In summary, ij}' — > p(770)7r is observed in ip' decays 
for the first time, and the branching fraction is measured 
to be Bii^' yo(770)7r) = (5.1 ± 0.7 ± 0.8) x IQ-^. A 
high mass excited p state at mass around 2.15 GeV/c^ 
is also observed with B{'4)' /9(2150)7r — + tt+tt^tt^) = 
(19.4±2.5t"i^) X 10"^. The results may shed light on the 
understanding of the longstanding "p7r puzzle" between 
J/V' and if)' hadronic decays. 
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